A precise determination of the Be mass from dynamical lattice QCD 
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We perform a precise calculation of the mass of the Be meson using unquenched configurations from the MILC 
collaboration, including 2+1 flavours of improved staggered quarks. Lattice NRQCD and the Fermilab formalism 
are used to describe the b and c quarks respectively. We find the mass of the Be meson to be 6.304(16)GeV. 
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1. Introduction 

One of the results anticipated from run II 
at the Tevatron is a measurement of the largely 
unexplored properties of the Bc{bc) meson. The 
current experimental determination comes from 
the CDF collaboration who quoted a value of 
6.4(4) GeV and from preliminary D9 data, 
where 231 events have been identified, giving the 
Be mass as Mb, = 5.95(37) GeV A suitably 
precise lattice calculation puts us in a position to 
make a prediction of this mass, rather than con- 
firm an experimental result. 

In the past, the two main theoretical methods 
for finding Mb, have been potential models |4l5j 
and lattice calculations in the quenched approxi- 
mation |6l7j . Both of these techniques have found 
agreement with the experimental values, but both 
have drawbacks; potential models inevitably de- 
pend on the form of the potential used, while the 
effect of quenching in a lattice calculation is esti- 
mated to be about 100 MeV. 

Now, however, it is possible to repeat the lat- 
tice calculation without the quenched approxima- 
tion. The key is to use the unquenched ensembles 
for lattice gauge fields from the MILC collabo- 
ration jQ. These have 2-1-1 fiavours of dynami- 
cal quarks. They use an improved gluon action. 
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and the Asqtad action for staggered light quarks, 
leaving discretisation errors of ©(a^a^). With the 
MILC ensembles, lattice calculations agree with 
experimental measurements for a wide range of 
hadronic quantities [5]. 



2. Method 

In bottomium (66) systems the typical velocity 
of the quarks is Vb ~ 0.1, while in charmonium 
(cc) Vc ~ 0.3. However, the unequal sharing of 
quark masses in the Be changes these values, giv- 
ing Vc ~ 0.5 and Vb ~ 0.04. With this in mind, 
we choose to use different formalisms to describe 
the 6 and c quarks. 

For the c quark we use the clover action, with 
the non-relativistic (Fermilab) interpretation jl()| . 
The clover coupling is adjusted to its tadpole im- 
proved tree-level value. For the 6 quark we use 
0{v'^) lattice NRQCD ^ with coefficients fixed 
at tree level. Our reasoning for these choices 
comes from the charmonium spectrum, repro- 
duced using Fermilab quarks ^2], and from the 
T spectrum, reproduced using the NRQCD ac- 
tion ^3]. Both of these calculations were per- 
formed on the MILC 2 -|- 1 fiavour ensembles 
which our calculations use. 

In both the NRQCD and Fermilab formalisms, 
the hadron energy at zero momentum is offset 
from the mass by an energy shift. In appropriate 
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mass differences, this shift cancels. Therefore, 



Mb, - i [M^ + Mr] 
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yield the differences in binding energies. (Here 
Mx and Ex denote masses and energies calcu- 
lated on the lattice respectively.) After comput- 
ing these differences on the lattice, we obtain our 
result for Mb^ by adding back [A4^ + M~f]/2 or 
Md^ + ■ These methods for extracting Mb^ 
will be referred to as the quarkonium baseline ^ 
and the heavy-light baseline method When 
results are quoted for the quarkonium baseline 
method, the "charmonium particle" used is the 
spin average of the rjc and J/ip- The absence of 
an experimental signal for rjb prevents us from 
doing this for the b quark as well. 



■light' u 1:2:3 ■ — |— 
'heavy' u 1:2:3 ,-Q- 
'uO' u 1:2:3 .-B- 
'fine.dat' u 1:2:3 : X 



® 



0.2 0.25 0.3 0.35 0.4 
sea quark mass ratio m|/mg 



Figure 1. Chiral limit of quarkonium and heavy- 
light baseline methods. Only statistical errors are 
shown. The m; /rUs ratio corresponds to the bare 
masses used for the dynamical quarks. 



3. Simulation Details 

In our calculation, the configurations included 
a single flavour at around the strange quark mass 
(ms) and two degenerate light flavours over a 
range of masses down to ms/5. The valence c and 
b quark masses were set using the kinetic masses 
of the Dg and T states. The lattice spacing is 
set on each ensemble through the 1S-2S radial 
splitting of the T system. But an important ad- 
vantage of calculations on the MILC ensembles is 
that any of several quantities shown in could 
have been used. 

In addition to the majority of our calculations 
performed at a « 0.12 fm we ran a single calcula- 
tion on a finer ensemble, at a w 0.09 fm. 

4. Results 

The result of calculating using different masses 
for the two (degenerate) light sea quarks is shown 
in Fig. 2] This plot also allows a comparison of 
the quarkonium and heavy-light baseline methods 
of Eqs. O and ©. 

We plot our result from calculating at the 
"fine" lattice spacing (a ^ 0.09/m). We find on 
these configurations that Mb^ = 6.309(3) GeV, 
(quoting only the statistical error). We also plot 



a single point where the tadpole improvement fac- 
tor (uq) for the b quark has been calculated using 
the mean link in Landau gauge as opposed to the 
fourth root of the plaquette (used for all other 
points). Changing only this quantity allows us to 
assess the infiuence of the relativistic corrections 
in the NRQCD action. We see that the effect is 
negligible, as expected. 

The effect of varying the valence quark masses 
on the kinetic masses of the J/ tp and T was used 
to set the error on the Be from these parameters. 
We find that this leads to a error of 10 MeV due 
to the uncertainty in the b quark mass and 5 MeV 
due to the uncertainty for the c. 

We also estimate a systematic error arising 
from higher order terms in the Fermilab quark ac- 
tion. For the heavy quarks, the NRQCD action 
includes the 0{v'^) terms at tree level, the Fer- 
milab quark action also includes them but with 
some mismatch. Because of the form of Eqs. ^ 
and |(5J we expect a more efficient cancellation 
of these differences using the quarkonium base- 
line method, where all quantities are in the heavy 
quark sector. Estimating the corrections, we get 
a ±10 MeV uncertainty for both methods, and an 
additional downward shift of 30 — 50 MeV for the 
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Figure 2. Mass of the Be from potential mod- 
els, Quenched lattice calculation and the present 
work (quarkonium baseline method). 



5. Conclusions 

This calculation shows once again the impor- 
tance of including dynamical quarks in lattice 
simulations. Having removed the ambiguity in 
setting the lattice spacing and fixing the quark 
masses, we are able to reach an unprecedented 
accuracy for a lattice calculation of the mass of 
the Be particle. This sets a well defined target 
for the current experimental studies. 
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heavy-light baseline method. Evidence for this 
shift can be seen in Fig. 

We give our final result using the quarkonium 
baseline method because of this shift. Our result 
is 6.304(3)(10)(11)(5)(3) GeV, where the first er- 
ror is statistical, the second is due to higher or- 
der corrections to the Fermilab quark action, the 
third is due to uncertainties in the valence quark 
masses, and the final number is an estimate of our 
chiral extrapolation error. Further error analysis 
is still underway, but combining the current errors 
in quadrature gives: 

Mb, = 6.304(16) GeV. (3) 

This result is shown in context in Fig. |5| We 
also show points due to Eichten and Quigg ^ and 
Kwong and Rosner who estimate results from 
a range of potentials. 

The result due to Shanahan [S] et al. is the 
previous lattice determination. Their calculation 
used the same method as this work, but was 
carried out in the quenched approximation. We 
achieve a significant reduction in the size of the 
systematic error from the removal of the quenched 
approximation. 
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